Bismuth qubits in silicon: the role of EPR "cancellation resonances" 
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We investigate theoretically and experimentally the electron paramagnetic resonance (EPR) spec- 
tra of bismuth doped silicon (Si:Bi) at intermediate magnetic fields, B f» 0.05 — 0.6 T. We identify 
a previously unexplored EPR regime of "cancellation-resonances"- where the non-isotropic part of 
ASzIz, the Ising part of the hyperfine coupling, is resonant with the external field- induced splitting. 
We show this regime has interesting and experimentally accessible consequences for spectroscopy 
and quantum information applications. These include reduction of decoherence, fast manipulation 
of the coupled nuclear-electron qubit system and line narrowing in the multi-qubit case. We test 
our theoretical analysis by comparing with experimental X-band (9.7 GHz) EPR spectra obtained 
in the intermediate field regime. 

PACS numbers: 03.67.Lx,03.67.-a,76.30-v,76.90.-|-d, 



Following Kane's suggestion [l[ for using phosphorus 
doped sihcon as a source of qubits for quantum comput- 
ing, there has been intense interest in such systems Q • The 
phosphorus system ('^^P) is appealing in its simphcity: it 
represents a simple electron-spin qubit S — ^ coupled to 
a nuclear-spin qubit I = via an isotropic hyperfine in- 
teraction AI.S of moderate strength (^ — 117.5 MHz). 
However, recent developments ;i o point to Si:Bi (bismuth 
doped silicon) as a very promising new alternative. Two 
recent studies measured spin-dephasing times of over 1 ms 
at 10 K which is longer than comparable (non-isotopically 
purified) materials, including Si:P 0, 0]. Another group 
implemented a scheme for rapid (on a timescale of ~ 100/xs) 
and efficient (of order 90%) hyperpolarization of Si:Bi into 
a single spin-state Q. 

Bismuth has an atypically large hyperfine constant ^ = 
1.4754 GHz and nuclear spin / = |. This makes its EPR 
spectra somewhat more complex than for phosphorus and 
there is strong mixing of the eigenstates for external field 
B < 0.6 T. Mixing of Si:P states was studied experimen- 
tally in [6\ , by means of electrically detected magnetic res- 
onance (EDMR), but at much lower fields B < 0.02 T. 
Residual mixing in SiiBi for _B = 2 — 6 T, where the eigen- 
states are > 99.9% pure uncoupled eigenstates of both 
and Sz, was also proposed as important for the hyperpo- 
larization mechanism of illuminated SiiBi 5]. In [4] it was 
found that even a 30% reduction in the effective param- 
agnetic ratio ^ (where / is the transition frequency) lead 
to a detectable reduction in decoherence rates. 

Below we present an analysis of EPR spectra for Si:Bi 
and test the results with experimental spectra. We iden- 
tify a series of regimes for which ^ — 0, explaining them 
in a unified manner as a series of EPR "cancellation reso- 
nances"; some are associated with avoided level-crossings 
while others, such as a maximum shown in ENDOR [7] 
spectra at i? w 0.37 T in 14] is of a quite different ori- 
gin. These cancellation resonances represent, to the best of 
our knowledge, an unexplored regime in EPR spectroscopy, 
arising in systems with exceptionally high A and /. They 



are somewhat reminiscent of the so-called "exact cancella- 
tion" regime, widely used in ESEEM spectroscopy 0, 13, 
but differ in essential ways: they affect both electronic and 
nuclear frequencies rather than only nuclear frequencies; 
they concern only the non-isotropic component of the in- 
teraction (and are thus not "exact"; indeed the B w 0.37 
T point is not even a full cancellation). They have impor- 
tant implications for the use of Si:Bi as a coupled electron- 
nuclear qubit pair: we show all potential spin operations 
may be carried out with fast EPR pulses (on nanosecond 
timescales) where in contrast, most operations for Si:P re- 
quire slower NMR pulses (on microsecond timescales). A 
striking spectral signature is reduced sensitivity to certain 
types of ensemble averaging, giving an analog to the ultra- 
narrow lines well-known in "exact cancellation" , as well as 
the reduction of decoherence. Further details are found in 

We model the Si:Bi spin system approximately by 
a Hamiltonian including an isotropic hyperfine coupling 
term: 



H = ijJqSz — ujqSIz + AS.i 



(1) 



where ujq represents the frequency of the external field and 
6 ^ uj /uo = 2.488 X 10^*^ represents the ratio of the nu- 
clear to electronic Zeeman frequencies. For / = |, 5 = ^ 
there are 20 eigenstates which can be superpositions of 
high-field eigenstates \ms,mi); but since [H,Sz + Iz] = 0, 
the \ms,mi) basis is at most mixed into a doublet |tos = 
±i, mj = TO =p i) with constant m = rUg + mj. One can 
thus write the Hamiltonian for each m sub-doublet as a 2 



dimensional matrix _ff,„ (where H„ 



- 



[■m+LJo{l+S)]c7z + {25-m^) 



-(-+2to<5c:;o)I (2) 



and where = ^ is the rescaled field, Oz , repre- 
sent Pauli matrices in the two-state basis |ms,TO/) = 
I ± ^,TO =F ^) and I is the identity operator. It be- 
comes clear that whenever to = — (I'o(l + 5), the quan- 
tum states become eigenstates of Gx- Thus at to w — lio; 



FIG. 1: The 20 spin energy-levels of Si:Bi may be labelled in 
alternative ways: (1) in order of increasing energy |2)...|20) 

(2) using the adiabatic basis |±,m) of doublets; levels belong- 
ing to the same doublet are shown in the same color (online). 

(3) by their asymptotic, high-field form \ms,mi). States |10) 
and |20) are not mixed. State |10) is of especial significance 
since it (rather than the ground state) is a favourable state to 
initialise the system in (experimental hyperpolarization studies 
Q concentrate the system in this state). Thus, in our coupled 
2-qubit scheme, state |10) corresponds to our (0)e(0)n state; in 
the same scheme, states |9) = (0)e(l)n and |11) = (l)e(0)n are 
related to 1 10) by a single qubit fiip (of either the electron or nu- 
clear qubits respectively) while for |12) = (l)e(l)n both qubits 
are fiipped. 
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the eigenstates, |±,m), assume Bell-like form: ( |±,m) 



^[|-i)|m+i)±| 



i)]). In contrast, "exact 
cancellation" results in a simple superposition of the nu- 
clear states, which also allows other types of manipulations 
0. Since > 0, only states with — 5 < m < can yield 
resonances where the field-splitting term \uJoaz is elimi- 
nated. In this case, they occur at ~ 0, A, ...4A, SA, cor- 
responding to applied field S = 0, 0.053, ...0.21, 0.26 Tesla; 
below we show that all ^ = points which are minima 
occur midway between these resonances. But the ^ = 
maximum at the wq = 7 A = 0.37 T resonance, and seen in 
experiments Q is shown to be of a different type. 

It is standard practice to represent two-state quantum 
systems using vectors on the Bloch sphere Q . We define a 
parameter i?^ = [m + wo(l + 5)] +25 — where Rm rep- 
resents the vector sum magnitude of spin x and z compo- 
nents. Denoting 9 as the inclination to the z-axis, cos 9m = 
[m + wo(l + 5)] /Rm and sin6'm = (25 — rr^Yl'^ jRm', then 
Eql5]can also be written: 



hra = Rm C0s9m^z + Rm Sa\9,n^x - 'j^i'^ + '^CjQmS)!. (3) 

Straightforward diagonalisation gives the pair of eigen- 
states, for each m, at arbitrary magnetic fields wq: 



|±,m) = a±| ± i m T + bfj T ^,to± ^) (4) 



FIG. 2: (a) Comparison between theory (Eq[6] and [7| (black 
dots) and experimental CW EPR signal (red online) at 9.7 GHz. 
Resonances without black dots above them are not due to Si:Bi; 
the large sharp resonance at 0.35 T is due to silicon dangling 
bonds while the remainder are due to defects in the sapphire 
ring used as a dielectric microwave resonator. The variation 
in relative intensities is mainly due to the mixing of states as 
in Eq[4l The variability is not too high but the calculated in- 
tensities are consistent with experiment and there is excellent 
agreement for the line positions, (b) Calculated EPR spec- 
tra (convolved with the 0.42mT measured linewidth); they are 



seen to line-up with the experimental spectra at / 



9.7 



GHz). The " cancellation resonances" are indicated by arrows 
and integers — m = 0, 1, 2, 3, 4.... The uo = 7^4 maximum of ^ 
indicated by the integer 7 corresponds to that shown in the < 2 
GHz ENDOR spectra of [J. 



where: 



sm ( 



and the corresponding eigenenergies: 



A 



1 



(1 + AiJomS) ± R„ 



(5) 



(6) 



In FiglHthe simple (but exact) expression EqE] reproduces 
the spin spectra investigated in e.g. [1 and fHS. Eqsl 
arc valid for all states except the unmixed m = ±5 states 
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(|10) and |20)). For \m\ = 5, there is no (Jx coupling: 
these two states are unmixed for all magnetic fields, thus 
a±5 = 1 and 6±5 = and Eq[S]siinphfies drasticahy to give 

For the doublets, the are the dominant coefficients 
at high-field. Then, the B — > oo limit corresponds to angle 
6m = 0, so — )■ 1 and b^^ — > and the states become 
uncoupled. The |m| < 4 cancellation resonances corre- 
spond to 9m — 7r/2 so — while 6^ = "'""Tl' "^^^ 
EPR emission transitions |-l-,r?T,) — | — ,^1 — 1) are dipole 
allowed at all fields. Their relative intensities 
/ (X 2/+--_i|(m, = ij^.lm; Since - 

i)P = i, variations in line intensities arise from mixing of 
the states. Thus, 

-^m^rn-l « = ^ (f COS 6*™) (f COS 6I,„_ i ) . 

If mixing is significant |+,rn) — > |-|-,m — f) transitions 
(of intensity /+) and \ — ,m — I) ^ | — ,™) transitions (of 
intensity /~), EPR-forbidden at high field, become strong, 
with relative intensities: 

/± = i(I±cos0„)(lTcos0™_i). (8) 

Forbidden lines disappear at high fields; as oo, 
one can see from Eq|5] that Ij-m^m-i ^ ^ — > since 
ItP cx at high fields. Eqs HI [7] and [8] are exact so are 
in complete agreement with numerical diagonalisation of 
the full Hamiltonian. 

In Fig|2^ we test the equations against experimental 
spectra at < 0.6 Tesla and microwave frequency 9.67849 
GHz. The long spin relaxation times at low temperatures 
means that the EPR spectra are easily saturated, compli- 
cating the analysis of the line intensities. We therefore used 
a temperature of 42 K so as to measure unsaturated res- 
onances. The shorter relaxation times at these elevated 
temperatures may be due to the presence of significant 
numbers of conduction electrons that are no longer bound 
to Si:Bi donors. We measure a very broad microwave ab- 
sorption centred on zero magnetic field (subtracted from 
Fig 12^) which we attribute to these conduction electrons. 
The comparison with experiment shows excellent agree- 
ment with the positions of the resonances, which are far 
from equally spaced in the low-field regime. Experimen- 
tal lines are found to be Gaussians of width 0.42mT; 
this was attributed in Q to the effects of ^^Si in this sam- 
ple; samples with enriched ^®Si are expected to give much 
narrower linewidths. 

For Figl2b, we generated the spectra, convolved with 
0.42mT Gaussians to obtain the EPR spectra of all lines 
(both allowed and forbidden) at all frequencies below 
f OGHz. We indicate the main dipole allowed lines as well 
as indicating the approximate position of the main reso- 
nances. The spectra show a striking landscape of tran- 
sitions which show maxima or minima where — 
and double-valued EPR resonant fields (i.e. transitions 
with EPR resonances at two different magnetic fields) . No 
Boltzmann factor has been included in the simulation. The 



nuclear field splittings are unresolved and extremely small; 
they do not affect line intensities significantly. To simplify 
our discussion, we neglect the tiny nuclear shifts cx luq6, 
but include them whenever spectroscopically significant. 

The well-studied 4-state S ^ 1/2, I = 1/2 Si:P sys- 
tem can be mapped onto a two-qubit basis. With a 20- 
eigenstate state-space, the Si:Bi spectrum is more complex, 
but we can identify a natural subset of 4 states, which rep- 
resents an effective 2-coupled-qubit analogue: 

112) ^(I).(I)„ 
Iff) ^(I),(0)„ 

|9) -^{0)e{l)n 
|fO) ^(0)e(0)„ 

As hyperpolarization initialises the spins in state |I0) and 
this state has both the electron and nuclear spins fully anti- 
aligned with the magnetic field, it can be identified with 
the (0)e(0)n state. The other states are related to it by 
either one or two qubit flips, just as in the Si:P basis. 

To have a universal set of gates for quantum informa- 
tion it is known to be sufRcient to be able to perform 
arbitrary single qubit manipulations and a control-NOT 
(CNOT) gate [ll]. In the two qubit system described here, 
arbitrary electronic qubit-only manipulations can be per- 
formed with radiation pulses exciting transitions between 
states 12 — 9 and 10 — ^ 11, while single nuclear-qubit rota- 
tions correspond to 12 -> 11 and 10 9. The CNOT gate 
(for example using the nuclear spin as a control qubit) is 
even simpler as it requires only a tt pulse connecting 12 — >■ 9 

The electronic flips are EPR allowed at all fields for both 
Bi and P donors so can be performed in a time on the order 
of 10 ns [3]. The nuclear transitions, however, require a 
slower, (of order microseconds) NMR pulse for Si:P. For 
Si:Bi, on the other hand, at the m = —4 resonance the 
nuclear and electronic transition strengths become exactly 
equal as may be verified by setting 6'_5 = and 6-4 = t^/2 
in Eqs.Q and Time-dependent calculations [12| show 
that the duration of a tt pulse is also equalized. 

The ujQ — A, 2A, 3A, AA resonances yield textbook level 
anti-crossings as well as "Bell-like" eigenstates. Were it to 
become possible to vary the external field sufficiently fast 
to produce sudden, rather than adiabatic evolution of the 
eigenstates it would be possible to transfer the former to 
the high-field regime. Unfortunately, ramping magnetic 
fields (up or down) sufficiently fast to violate adiabaticity, 
though not impossible, would require some of the fastest 
magnetic field pulses ever produced (eg lO^T/s obtained 
by [lO|). However, we show that adiabatic magnetic 
field sweeps already achievable by ordinary laboratory 
magnetic pulses (1 — 10 T/ms) suffice to already achieve 
new possibilities. 

The frequency minima at 5-8 GHz: fields at which 
df/dB = are expected to lead to a reductions of deco- 
herence, since sensitivity to magnetic fluctuations is min- 
imised; a measurable reduction has been seen by varying 
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FIG. 3: The simulated EPR spectra for pairs of interacting 
Si:Bi donors shows ultra-narrow lines at cancellation resonances, 
somewhat analogous to "exact cancellation" in ESEEM spec- 
troscopy. The electronic spins are coupled by an exchange term 
JS1.S2. The left panel shows the effect of a typical J — 0.5 
GHz coupling. A splitting (of order J) appears in general, but 
near the m = —4, —3 resonances, this is suppressed. The graphs 
on the right show the calculated signal from a sample of spins 
with a distribution of J, where average (J) = 0.3 GHz, variance 
aj — 0.3. At the cancellation resonances, despite the ensemble 
averaging, the lines remain narrow. 



the ratio df /dB by 30 — 50%. In Fig.® we see several tran- 
sitions have a minimum frequency. These minima (in effect 
of Rjn + Rm-i) occur for: 



cos Or, 



■ cos ^ra-\- 



(9) 



Thus Qm = — Ora~\^ the conscqucnce is that the min- 
ima He exactly midway in angular coordinates between 
cancellation resonances. For example, for the 12 ^ 9 
hne (| + ,m = —3) | — ,m = —4)) the minimum is at 
Wo = 2|i ^ 3.57yl so B ~ 0.188 T. Here, the m = -4 dou- 
blet has passed its resonance point at 0.21 T (for which 
9-4 — n/2) by an angle 4> = arccos and the m = — 3 
resonance is at an equal angular distance before its reso- 
nance at ~ 0.16 T: thus 6'_4 = 7r/2-h(;!) while 6-3 = iT/2-(j). 
Both doublets are quite close to the Bell-like form. 

Line narrowing: an interesting and unexpected conse- 
quence has applications to studies with larger numbers of 
spins. A pair of Bi atoms, interacting via a spin-exchange 
term of the form JS1.S2 will result in splitting of the EPR 
spectral lines (with an energy splitting of order J) . How- 
ever this is suppressed near the cancellation resonances. 
Analogously to "exact cancellation", this makes the sys- 
tem less sensitive to ensemble averaging. For exact can- 
cellation, this means the averaging over different orienta- 
tions in powder spectra; here it means magnetic pertur- 
bations including spin-spin interactions. Fig|3] (left panel) 
plots the signal for J = 0.5 GHz for a single pair of Si:Bi 
atoms and clearly shows the line splitting away from the 
resonances. The right panels show the effects of averaging 
many spectra each corresponding to different J (with an 
average (J) = 0.3 GHz and width aj = 0.3 GHz). While 
typical spectra show a broad feature of width ^ ctj, at 
the cancellation resonance the line width remains strikingly 
narrow (close to the single atom line width). 

A frequency maximum at wq = 7 and B ~ 0.37 T is 
marked with a 7 in FigI5]D. We can show that a df /dB = 



point which is a maximum (in effect of Rm — Rm-i) implies: 

COS^m = C0s6'm„i. (10) 

This condition does not correspond to the elimination of 
the field splitting terms; instead it implies 0_3 = 9-^ = 
7r/4, thus equalizing the Bloch angle for the associated en- 
ergy levels. In this sense it is somewhat different to the 
other cancellation resonances; nevertheless, it still provides 
a df/dB = point and thus so me p otential for reducing 



broadening and decoherence. In [1^, it is shown that the 
luq — 7 resonance offers new possibilities for copying and 
storing qubit states. At = 5, the most drastic cancella- 
tion resonance occurs, since both cr^ and terms in Eq|2] 
are eliminated, leaving only the isotropic term. Although 
there is no df/dB = or line narrowing here, there is a pos- 
sibility of driving, by a second order process, simultaneous 
qubit rotations e.g. |0)e|l)„ — )■ |l)e|0)„ see |12[. 

Conclusions: In the intermediate-field regime {B ~ 
0.05 — 0.6 T) the exceptionally large values of A and / for 
Si:Bi generate a series of cancellation resonances. They are 
associated not only with level crossing structures but also 
with more subtle and not previously studied effects: both 
line broadening and decoherence effects may be reduced; 
also, if the electronic and nuclear spins of Si:Bi are used as 
a 2-coupled qubit system, the cancellation resonances allow 
a universal set of quantum gates to be performed with fast 
EPR microwave pulses, eliminating the need for slower ra- 
dio frequency addressing of the nuclear qubit. One scheme 
would envisage the following stages: (1) hyperpolarization 
of the sample into state |10) (in which the 2-qubit sys- 
tem is initialized as |0)e|0)„) at i? « 5T. (2) A magnetic 
field pulse (~ 10 T/ms, of duration lower than decoherence 
times) would reduce _B to ~ O.IT. (3) As the pulse ramps 
up, a series of EPR pulses would execute a series of gates 
and operations on the system. (4) As the magnetic pulse 
decays, the system is restored to the high-_B limit, leaving 
it in the desired superposition of |0)e|0)„, |l)e|l)n, |0)e|l)n 
and |l)e|0)„ basis states. Thus, given the capability to 
rapidly (< 1ms) switch from the high to intermediate field 
regime, Si:Bi confers significant additional possibilities for 
quantum information processing relative to Si:P. 
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